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■ Abstract. We argue that an expanding quark-gluon plasma has an anomalous 
lO ! viscosity, which arises from interactions with dynamically generated colour fields. The 

anomalous viscosity dominates over the collisional viscosity for large velocity gradients 
or weak coupling. This effect may provide an explanation for the apparent near perfect 
— I liquidity of the matter produced in nuclear collisions at RHIC without the assumption 

■ that it is a strongly coupled state. 
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O ■ 1. Introduction 



Measurements of the anisotropic collective flow of hadrons emitted in noncentral 
collisions of heavy nuclei at the Relativistic Heavy Ion Collider (RHIC) are in remarkably 
^ ■ good agreement with the predictions of ideal relativistic fluid dynamics. The comparison 
between data and calculations indicates that rj cannot be much larger than the 
^ ■ postulated lower bound r/min = s/An [U |2[ [3]. This have led to speculations that 
. the matter produced at RHIC is a strongly coupled quark-gluon plasma (sQGP). The 
possible microscopic structure of such a state is not well understood at present B]. 

Here we present an alternative mechanism that may be responsible for a small 
viscosity of a weakly coupled, but expanding quark-gluon plasma. The new mechanism 
is based on the theory of particle transport in turbulent plasmas [6l [7] . Such plasmas are 
characterized by strongly excited random field modes in certain regimes of instability, 
which coherently scatter the charged particles and thus reduce the rate of momentum 
transport. The scattering by turbulent fields in electromagnetic plasmas is known to 
greatly reduce the viscosity [8l|9] of the plasma. Following Abe and Niu [Qj, we call the 
contribution from turbulent fields to transport coefficients "anomalous" . 

The sufficient condition for the spontaneous formation of turbulent, partially 
coherent fields is the presence of instabilities in the gauge field due to the presence of 
charged particles. This condition is met in electromagnetic plasmas with an anisotropic 
momentum distribution of the charged particles p^, and it is known to be satisfied 
in quark-gluon plasmas with an anisotropic momentum distribution of thermal partons 
[m [121 [13]. 
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The turbulent plasma fields induce an additional, anomalous contribution to the 
viscosity, which we denote as rjx- This anomalous viscosity decreases with increasing 
strength of the turbulent fields. Since the amplitude of the turbulent fields grows with 
the magnitude of the momentum anisotropy, a large anisotropy will lead to a small 
value of Tjx. Because the relaxation rates due to different processes are additive, the 
total viscosity is given by ri~^ = rfj^ + r^^^, where r/c is the collisional shear viscosity. 
This implies that r/A dominates the total shear viscosity, if it is smaller than r/c- In that 
limit, the anomalous mechanism exhibits a stable equilibrium in which the viscosity 
regulates itself: The anisotropy grows with 77, but an increased anisotropy tends to 
suppress t^a and thus r/ ~ r/A- 



2. Anomalous viscosity 

Here we give a heuristic derivation of ?7a. For more systematic derivations, see references 
[m [15] . According to classical transport theory, the shear viscosity is given by, 

7] ^ ^npAf, (1) 

where n denotes the particle density, p is the thermal momentum, and Af is the mean 
free path. For a weakly coupled quark-gluon plasma, n 5T^ and p ~ 3T. The 
mean free path depends on the mechanism under consideration, r^c is obtained by 
expressing the mean free path in terms of the transport cross section A^^"* = (natr)"^. 
Using the perturbative QCD expression for the transport cross section in a quark-gluon 
plasma yields the result, which agrees parametrically with the result for the collisional 
shear viscosity in leading logarithmic approximation [16]. The anomalous viscosity is 
determined by the same relation ([ID for r], but the mean free path is now obtained by 
counting the number of colour field domains a thermal parton has to traverse in order 
to "forget" its original direction of motion. If we denote the field strength generically by 
B"" (a denotes the colour index), a single coherent domain of size causes a momentum 
deflection of the order of ~ gQ°'B°'r^, where Q° is the colour charge of the parton. 
If different field domains are uncorrelated, the mean free path due to the action of the 
turbulent fields is given by 

Ar> = ,.„((p7Ap)^)~^3^^. (2) 

The anomalous shear viscosity thus takes the form: 

np^ 9sT^ 

The argument now comes down to an estimate of the average field intensity {B^) 
and size rm of a domain. We first note that the size is given by the characteristic wave 
length of the unstable field modes. Near thermal equilibrium, the parameter describing 
the influence of hard thermal partons on the soft colour field modes is the colour-electric 
screening mass ~ gT. Introducing a dimensionless parameter ^ for the magnitude 
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of the momentum space anisotropy [13] , the wave vector domain of unstable modes is 
fc^ < ^m"^. Thus Tin ~ ^"^/^((^T)"^. The exponential growth of the unstable soft field 
modes is saturated, when the nonlinearities in the Yang-Mills equation become of the 
same order as the gradient term: g\A\ ~ k, which implies that the field energy in the 
unstable mode is of the order of g'^{B'^) ~ fc"^ ~ ^^^d- The denominator in ([3]) thus has 
the characteristic size, at saturation: 

g'Q^B')r^r^e"<-e'\9Tf. (4) 

If our analysis, which requires confirmation by numerical simulations, is correct, it gives 
the following relation for the anomalous viscosity: 

~ • (5) 

We conclude that r^A will be smaller than ?7c, if the coupling constant g is sufficiently 
small and the anisotropy parameter ^ is sufficiently large. 

Reference [I3] uses the following parametrization of the anisotropic momentum 
distribution: 



\2 



/(P) = /o [^Jp' + e(p • ny) ^ foip) - -f^M^ ± /o) , (6) 

where /o is the thermal equilibrium distribution. For the special case of boost-invariant 
longitudinal flow, we obtain a relation between rj and ^, which takes the form (for a 
massless parton gas): 



s T 
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(7) 



Combining ([5]) and and setting t] = rjA, we obtain: 
s \g^\\/u\J 

The viscous contribution to the stress tensor is proportional to r^lVw]. For the 
collisional viscosity, this implies that the stress tensor grows linearly with |Vm|. The 
anomalous shear viscosity ([H]), on the other hand, is a decreasing function of the 
velocity gradient; its contribution to the stress tensor grows like IVmP*^^ for our scaling 
assumptions. The unusual dependence of r]A on |Vm| certainly justifies the term 
"anomalous viscosity" . 

The different dependence of the collisional and the anomalous viscous stress tensors 
on the velocity gradient is shown schematically in figure [H For very small gradients 
the linear dependence of the collisional viscous stress tensor dominates, but for larger 
velocity gradients the lower power associated with the anomalous shear viscosity will 
assert its dominance. The precise location of the cross-over between the two domains 
depends on the value of the numerical constant cq, but we can deduce from ([H]) that 
the cross-over point shifts to lower values of |Vm| with decreasing coupling constant g. 
We also show in the figure the effect of choosing a different power n in the scaling law, 
g^iB"^) = bog^T^^"^ for the energy density of the turbulent colour fields. Finally, we 
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note that the Kovtun-Son-Starinets (KSS) bound 77 > s/47r [3] does not apply to the 
anomalous viscosity. 

We have discussed the anomalous viscosity in an anisotropically expanding quark- 
gluon plasma. By reducing the shear viscosity of a weakly coupled, but expanding 
quark-gluon plasma, this mechanism could possibly explain the observations of the RHIC 
experiments without the assumption of a strongly coupled plasma state. 




Vu/T 



Figure 1. Schematic representation of the dependence of the coUisional and 
anomalous viscous stress on the velocity gradient. The coUisional viscous stress is 
shown by the linearly rising line; the anomalous viscous stress is shown by the curved 
lines for three scaling exponents of the turbulent colour field energy (n = 1.5,2,2.5). 
The solid lines indicate the dominant source of viscous stress in different regions of the 
scaled velocity gradient |Vu|/T. 
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